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In this Letter we map out the mean field energy potential landscape of fermion pairing states 
with population imbalance near broad Feshbach Resonances. We apply the landscape to investigate 
the nature of phase separation, when the Hilbert space is subject to the constraint of constant 
population imbalance. We calculate the scattering length dependence of the critical population 
imbalance for various phase separated states across Feshbach resonances. 



Recently, cold atom pairing states with population dif- 
ferences were studied experimentally by the MIT and 
the Rice University cold atom groups 0, 0- Fermion 
pairing with population imbalance has been a fascinat- 
ing subject for a longtime and was studied in vari- 
ous contexts 0, 0, 01 M 01 13- For cold atoms, pair- 
ing occurs between two hyperhne states participating 
in Feshbach resonances, as indicated by experiments a 
few years agojjj. Cold atom superfluids across Feshbach 
resonances with population imbalance have also been 
studied theoretically, to some extent; different phases 
were proposed, based on a variety of arguments or 
calculations [Tol HH H^ |. Most recently, cold atom su- 
perfluids with population imbalance in hnite traps were 
addressed 0. It is generally believed that such super- 
fluids exhibit different behaviors on cither side of reso- 
nances. In the BCS limit, where the scattering length 
is negative and small in magnitude, the superfluid un- 
dergoes a first-order phase transition to a LOFF state 
(pairing between nonzero total momentum), followed by 
a second order phase transition to a partially polarized 
normal state, as the difference in chemical potentials be- 
tween the two spin species is increased |5|,|6j. On the other 
hand, on the far BEC side where the scattering length is 
positive and small, one has a continuous transition to a 
superfluid state coexisting with a Fermi sea composed of 
the excess spin species. 

One of the unique features of cold atoms is that re- 
laxation processes including the (global) population im- 
balance relaxation are very slow compared to the exper- 
imental measurement time period. Therefore, for cold 
atoms the population imbalance Pi m can be considered 
to be independent of /, the energy splitting between two 
fermions with the same momentum. This is opposed 
to traditional solid state, many-body systems where the 
chemical potential difference between the two species is 
a unique function of I, the Zceman splitting energy, in 
equilibrium. For cold atoms, the relevant ground state 
is therefore the state of lowest energy in the constrained 
Hilbert space defined by fixed total particle number N 
and population imbalance Pi m , not the lowest energy 
states in the whole Hilbert space for a given total den- 



sity. Motivated by the above observation and the recent 
experiments on cold atom pairing, in this article we map 
out the energy potential^ V) landscape of all lowest en- 
ergy homogeneous states in different subspaccs specified 
by Pim- (We leave out discussions on LOFF or other 
more exotic competing translational-symmetry-breaking 
states in this article.) The landscape not only provides 
important information on the energetics of various states 
but is also critically valuable for the understanding of 
phase separation phenomena (see below). Furthermore, 
it paves the way for future studies of the dynamics of 
superfluids in the limit of slow relaxation. As the first 
application, we study the phase separation with the help 
of the landscape and construct phase separated states at 
a given population imbalance Pi m (and with a given den- 
sity p). Finally we investigate critical population imbal- 
ance for various phase separated states across Feshbach 
resonances. 

The phase separation phenomenon discussed here only 
occurs when we consider states in a Hilbert subspace with 
a conserved population imbalance. These phase sepa- 
rated states generally do not correspond to ground states 
in the whole Hilbert space for a fixed total number of 
particles. In a population-imbalance-conserved subspace 
phase separation occurs at any energy splitting I. The 
main results are summarized in Fig.l. Unlike in usual 
superconductors where a state prepared at an arbitrary 
Pi m and I generally is driven toward an equilibrium state 
with Pi m — Pf m (I), a cold atom state can only settle 
with the lowest energy state in the Hilbert subspace of 
a given initial population imbalance Pj m . Phases in the 
inset of Fig. ^ ( a ) an d (b) are for those lowest energy 
states with a given population difference Pj m and den- 
sity p in a given field I. In Fig. ^ ( a ) an d (b) we present 
the scattering length dependence of critical population 
imbalance P c i, P C 2. 

The model we employ to study this subject is the 
standard one-channel model which is suitable for the 
discussion of superfluids near broad Feshbach reso- 
nances. The homogeneous states under consideration 
are BCS states with excess quasi-particles accommodat- 
ing unpaired atoms or the population imbalance, i. e., 
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FIG. 1: The critical population imbalance as a function of 
scattering length. In (a)[(b)], P c \(Pc2) versus (a s fc/) _1 for 
negative(positive) scattering length a s . P c i{P C 2) is the critical 
value below(beyond) which phase separation takes place(see 
also inset). Here N stands for a normal state, PS for a phase 
separation state, SF™ for a superfluid with a finite uniform 
population imbalance, and FP for a fully polarized state. In 
(b), the phase separated state is constructed using a super- 
fluid with a finite population difference and a fully polarized 
state. In inset, we show phases on the Pi m — I plane. Shaded 
areas are for phase separated states; the dashed lines are the 
schematic equilibrium curves P; m = P^ (I) for a given chem- 
ical potential, towards which the system is driven if an ex- 
change of particles with a reservoir is allowed. 



riker t^CS) (7^. ^ is the creation operator of quasi- 
particles). These states are analogues of the breached- 
pair states discussed previously^. Notice here \BCS) is 
a BCS state with a P im -dependent gap A(P im ) and 7^ 

is the quasiparticle creation operator. The corresponding 
gap is self-consistently determined by the following gap 
equation in the presence of population imbalance 
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-k = ^ / 2m ' ^k = e k - Ek — y " 
njj 1 is the occupation number of quasi particles at state 
k, spin-up and Xr is related to a s by XrV = —4ira s /m. 
Furthermore, only the lowest energy quasi-particle states 
are occupied so as to minimize the kinetic energy for 
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FIG. 2: The momentum distribution function f(k) of spin- 
up(solid line) and spin-down atoms(dashed line) in homoge- 
neous states studied in this article, (a) For negative scat- 
tering length and positive chemical potential; (b) for posi- 
tive scattering length and negative chemical potential. Inset: 
schematic of the occupation of quasi-particles. 
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which is unity only when \ey. — fi\ is less than the cut-off 
energy <5£. The energy of quasi-particles is Ek — I and can 
be either positive or negative; the distribution function 
here is not the usual equilibrium one. The total number 
of these quasiparticles precisely yields the population im- 
balance as indicated in Fig. El 
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Each state under consideration here is therefore the low- 
est energy homogeneous state in the corresponding con- 
strained subspace. Finally, the total number of particles 
is given as 
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The mean-field £V , ft —< Hbcs > —fJ-N (or the free 
energy at T = but with Pi m fixed) is 



^E^k-'H.t + E^k-^k 

k k 



A\ 
2s k ' 



A 2 



(5) 



Taking into account Eq. @ one can also easily obtain the 
energies of states with given N and Pi, n . 

The EV landscape of homogeneous states at different 
scattering length has been obtained. When the scatter- 
ing length is negative and the magnitude is small, we 
find a set of homogeneous solutions A(Pj m ) to the gap 
equation for a given population imbalance Pi m . The self- 
consistent BCS gap decreases as the magnetization Pi m 
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increases and becomes zero at a critical value as seen in 
Fig. 03 (b) . The critical value of Pi m is independent of the 
energy splitting as also indicated by vertical lines in Fig. 
13 (a), (c) and (e). The cold atoms in this limit would 
have a continuous phase transition to a normal state if 
the spatial phase-separation were prohibited. The EV of 
states with a given chemical potential and in the pres- 
ence of a given energy splitting is plotted as a function 
of the population imbalance Pj TO in Fig. |21 Firstly, the 
BCS state represented by the solution at Pi m = point 
becomes degenerate with a normal state (a solution with 
A = 0) at a critical magnetic field I c — Ao/v2- Sec- 
ondly, a local minimum (a metastable normal state) and 
a maximum (the Sarma solution) appear when the en- 
ergy splitting is between / = Aq/2 and I = Ao; here 
the normal and BCS state are interpolated by the Sarma 
solution as expected in the conventional BCS theory Q. 

To find the ground state for a given conserved popula- 
tion imbalance, we examine the EV landscape of the ho- 
mogeneous states. We find that at any given energy split- 
ting the energy curve is concave (down) when the pop- 
ulation imbalance is small. For a given population dif- 
ference, one constructs a mixed state which involves the 
BCS state and a normal state located at the end point of 
the concave part of the curve. These two states involved 
in the construction correspond to two degenerate states 
at the first order phase transition point and always have 
the same chemical potentials for each fermion species. 
For a given density p and Pi m , X the fraction of the 
BCS state and fi should be determined self-consistently 
by the following conditions, 
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P = Pbcs((*)X + Pn{p, 5p = 4(m))(1 ~ x ), 

P lm = Pun(p, Sp = J c (/i))(l - X). (6) 

Here Pbcs(N) is the density of a BCS (normal) state 
when the chemical potential difference between two 
species Sfj, is equal to the critical energy splitting I c ; Pi m 
is the corresponding population imbalance of the normal 
state. The phase separation occurs when the population 
imbalance is smaller than a critical one. The critical value 
corresponds to the population imbalance at the end point 
of the concave-downwards part of the curve, or the solu- 
tion to the above equation at X — 0. Above that critical 
value the homogeneous states regain stability. By car- 
rying out similar analysis at different scattering length, 
we also obtain the scattering length dependence of the 
critical population imbalance (Fig. ^a)). 

When the scattering length becomes positive and es- 
pecially when the chemical potential becomes negative, 
the EV landscape experiences an important qualitative 
change. The part of curve for the small population imbal- 
ance becomes convex (or concave up) instead of concave 
(down). On the other hand, beyond a critical imbalance 
population the curve again is concave. A transition from 



FIG. 3: In (a), (c) and (d) we show the energy potential as a 
function of the population imbalance Pi m at different energy 
splitting / but with a fixed chemical potential, (a) is for a 
negative scattering length (a a fc^) _1 = —1.40 and /1 = 0.98e%; 
(c) for positive scattering length near resonance (ajij.)" 1 = 
1.27 and fi — — 1.68e^; (d) for a small positive scattering 
length (a s fcf ) _1 = 1.56 and u = — 2.QGe° F . k° F is defined as 
the Fermi momentum for the density No /V at zero population 
imbalance Pi m — and e F is the corresponding Fermi energy. 
In (b), we show the variation of A and iV as a function of Pi m 
in the homogeneous solutions to Eq. (l)-(4) (with the same 
scattering length and chemical potential as in (a)). Notice in 
(c) and (d), all the states to the right of the vertical dashed 
line correspond to fully polarized states. Here N and Pi m are 
expressed in unit of A^o and is in unit of 0.75 2 ^ 3 e° F . 

convex to concave as the population imbalance increases 
is a general feature of the landscape when the chemical 
potential is negative and again is independent of the en- 
ergy splitting. 

As the energy splitting increases, the global minimum 
along the curve shifts continuously from Pi m = point 
into a uniform superfluid with a finite population imbal- 
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ance. This signifies a continuous transition from a BCS 
state to a uniform superfluid with population imbalance 
when the chemical potential becomes negative. Further 
increasing the energy splitting leads to an additional lo- 
cal minimum and local maximum (see Fig. yj (c)). The 
new local minima located by the vertical dashed line rep- 
resent fully polarized states and therefore the pairing am- 
plitude associated with these states is zero, i.e. A = 0. 
At an upper critical value, we obtain two degenerate so- 
lutions located at two ends of the concave part of the 
curve, one being a homogeneous superfluid with popu- 
lation imbalance and the other a fully polarized state. 
One uses these two states to construct a phase separated 
state since again the chemical potential for each fermion 
species can be shown to be the same. Phase separation 
occurs when the population imbalance exceeds a criti- 
cal value specified by the left end of the concave curve 
in Fig[!Jc) (i.e. the superfluid with uniform population 
imbalance). Note that in this case the density variation 
along a given curve is substantial. At a smaller posi- 
tive scattering length or deep into the BEC side of the 
resonance, we find the concave-down section of a curve 
moves toward higher population imbalance and finally 
disappears in the EV landscape (shown in Fig. 01(d)). 
Homogeneous superfluids with population imbalance are 
always stable in this limit. 

Before concluding we would like to make two com- 
ments. The first one is about the role of population 
imbalance relaxation (PIR). In the presence of PIR, the 
population imbalance takes a unique equilibrium value 
P ( ^(/); phase separation of cold atoms in this limit (with 
the total number density fixed) needs a fine tuning in the 
energy splitting /. If the difference between the density 
of a BCS state and that of a normal state at the same 
chemical potential is negligible, the phase separation only 
occurs when J takes a critical value. Indeed, as illustrated 
before 12] the phase separation takes place in the vicinity 
of a critical splitting and the interval becomes to be ex- 
ponentially small in the weakly interacting limit. In the 
view of this, the phase separation observed in Ref.Q, Q 
in the extreme quantum limit (i.e. extremely slow imbal- 
ance relaxation) is distinct from phase separation phe- 
nomena in the opposite limit where only equilibrium 
states are considered. Meanwhile, in general the time- 
dependence of population imbalance can be expressed as 
Pim(t) = P^(0)exp(-t/T)+P^(I)(l-exp(-t/r)) (here 
r is the imbalance relaxation time.). The temporal be- 
havior of phase separation therefore does depend on the 
energy splitting /, though practically such dependence 
is always suppressed if the measurement time is much 
shorter than the relaxation time r. 

The second remark is on the effect of confinement or 
finite size. In the presence of smooth potentials, the spa- 
tial variation of particle density effectively provides an 
ensemble of cold atoms superfluids at different chemi- 
cal potentials. The phase separation occurs in a region 



where the self-consistent chemical potential difference S/j, 
(determined by the population imbalance) is equal to the 
critical value 8[i c for the local self-consistent chemical 
potential (i(r). However, if the confinement potential 
is hard-wall like, then the chemical potential difference 
could be self- consistently pinned in the vicinity of the 
critical value 5[i c for the bulk density. 

In conclusion, we have mapped out the £ V landscape 
for fermion states with population imbalance. This re- 
sult might be the blueprint for the future studies of the 
dynamics of various competing states near Feshbach reso- 
nances. For instance it can be used to understand the en- 
ergy released, during the phase separation, from a homo- 
geneous state initially prepared and the evolution of ini- 
tial states. The phase separation of fermion pairing states 
in a population imbalance conserved subspace turns out 
to be robust and needs no fine tuning in the energy split- 
ting. This is also the limit where experiments were car- 
ried out and we believe the results obtained here will help 
to identify new interesting experimental opportunities in 
the future. We would like to thank Mike Forbes, Randy 
Hulet, Wolfgang Ketterle, Tony Leggett, and Boris Spi- 
vak for numerous discussions. This work is supported by 
a grant from the office of the Dean of Science at Uni- 
versity of British Columbia and a Discovery grant from 
NSERC, Canada. FZ is an A. P. Sloan fellow. 
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